On the origin of solar wind. Alfven waves induced jump of coronal temperature 



T.M. MishonovQ M.V. StoevQ and Y.G. ManevfS 

Department of Theoretical Physics, Faculty of Physics, 
University of Sofia St. Clement of Ohrid, 
5 J. Bourchier Boulevard, Bg-1164 Sofia, Bulgaria 
(Dated: February 3, 2008) 

Absorbtion of Alfven waves is considered to be the main mechanism of heating in the solar corona. 
It is concluded that the sharp increase of the plasma temperature by two orders of magnitude is 
related to a self-induced opacity with respect to Alfven waves. The maximal frequency for prop- 
agation of Alfven waves is determined by the strongly temperature dependent kinematic viscosity. 
In such a way the temperature jump is due to absorption of high frequency Alfven waves in a nar- 
row layer above the solar surface. It is calculated that the power per unit area dissipated in this 
layer due to damping of Alfven waves blows up the plasma and gives birth to the solar wind. A 
model short wave-length (WKB) evaluation takes into account the l// 2 frequency dependance of 
the transversal magnetic field and velocity spectral densities. Such spectral densities agree with 
old magnetometric data taken by Voyager 1 and recent theoretical calculations in the framework 
of Langevin-Burgers MHD. The present theory predicts existence of intensive high frequency MHD 
Alfven waves in the cold layer beneath the corona. It is briefly discussed how this statement can 
be checked experimentally. It is demonstrated that the magnitude of the Alfven waves generating 
random noise and the solar wind velocity can be expressed only in terms of satellite experimental 
data. It is advocated that investigation of properties of the solar surface as a random driver by 
optical methods is an important task for future solar physics. Jets of accretion disks are speculated 
as a special case of the wind from magnetized turbulent plasma. 

PACS numbers: 52.35.Bj, 52.50.Sw, 96.60.P-, 94.30.cq 



I. INTRODUCTION 

The coronal heating mechanism^ is one of the most 
perplexing longstanding unresolved problems of contem- 
porary physics^ This is the reason why even some con- 
densed matter theorists are involved in the revealing 
of this mystery. The lack of experimental data makes 
distinguishing between the various theoretical models 
difficult. In such a situation it is possible to apply 
purely aesthetic criteria for a natural description of some 
well-known facts. One such fact is that the increase 
of plasma temperature in the transition region by two 
orders of magnitude^ occurs very rapidly - there is a 
smeared jump, which serves as a starting point for cre- 
ation of the solar wind. We wish to emphasize that 
complete magneto-hydrodynamic (MHD) simulations of 
Alfven waves generated by a random driver also give 
a sudden increase of the temperature as a function of 
height £ The purpose of the current work is to give a qual- 
itative explanation of the observed temperature jump 
and to perform a model evaluation of its order of mag- 
nitude. Finally, we will analyze some possible future ex- 
perimental observations described by the presented the- 
oretical scenario. For example, the model predicts exis- 
tence of intensive high frequency Alfven waves beneath 
the coronal temperature jump. 

II. SCENARIO 

In order to make a proper assessment for the magnetic 
field spectrum comparable with magnetometric data 



taken by Voyager 1— we need to switch to a frequency- 
dependent regime and perform a time averaging of the 
Fourier transformed wave component B(t) w B x x + B y y 
perpendicular to the constant magnetic field Bo = Bq^ z z 
along which the z-axis is chosen 
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B(t) = 2jexp(-iw„t)B aJrl , uj n = —n, (1) 
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B w „ = — / exp(iu> n t)B(t) dt, n = 0, ±1, ±2, . . . 
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As our purpose is to present a model evaluation, from 
now on we will consider only one of the fluctuating wave 
components, so the plane magnetic field indices will be 
omitted. Fourier analysis of the signal accumulated for 
a time interval At gives us the opportunity to observe 
the magnetic field spectral density Q(t), which is tightly 
bound to the time averaged square of the wave magnetic 
component 
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The thus-defined spectral density 

I = Lj/2n, (/ 2 - f{) At » 1 (3) 

determines the magnetic field energy density w and the 
Alfven waves energy flux S as a function of the circular 
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In interplanetary space the spectral density $ is mea- 
sured in typical units nT 2 /Hz. Alfven waves energy den- 
sity depends on the integrated spectral density and the 
Alfven speed Va- In general the integration should be 
taken over all possible wave frequencies / £ (0, oo) un- 
less the given specific power spectrum imposes an intro- 
duction of a low or high frequency cut-off. In the sug- 
gested model we use the natural cut-off frequency set by 
the condition for waves existence loata = 1, when the 
Alfven waves frequency wa exactly equals the attenua- 
tion coefficient ta _1 . Therefore in the investigation that 
follows we exclude the influence of any extremely high- 
frequency Alfven waves that would be absorbed by the 
medium immediately after their generation under the so- 
lar surface. Hence for heating in the corona we only take 
into account those waves, whose characteristic life time 
exceeds the inverse value of the cut-off frequency u> c 
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Molecular kinetic theory! determines the temperature- 
dependence of the cut-off frequency, being a function 
of the kinematic viscosity v. A more detailed analysis 
shows that the sum of kinematic and magnetic viscosity 
is relevant in the damping of the Alfven waves^, but 
the magnetic viscosity is much smaller and negligible for 
sufficiently high plasma temperatures. This is because 
the kinamatic viscosity is cx T 5 / 2 while the resistivity is 
T- 3 / 2 . For more detailed discussion see, for example, 
the work by Erdelyi and Goossensi^and Ref. [111. For low 
density coronal plasmas, electromagnetic emission and 
thermal conduction are negligible for the energy balance 
close to the temperature jump. In our model we neglect 
the difference between the electron and ion temperatures 
which are of the same order T e « T p . However we wish 
to emphasize that heating is due to ion viscosity and our 
consideration gives a natural explanation why the pro- 
ton temperature is higher than electron one T p > T e . 
It is because viscous friction heats the ions. Since we 
pursue just a qualitative estimation, with a logarithmic 
accuracy in the final results for the solar wind velocity 
presented below, we can neglect the slight temperature 
variations in the protons' Coulomb logarithm L p as well 
as the influence of the electron temperature dependence 
on the Debye radius a. Precise calculations including 
such a dependence would only lead to minor corrections 
that would not change the order of magnitude of the final 
outcome 
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According to valuable analysis on the experimental data 
obtained by the Voyager 1 magnetometer—, the magnetic 
field spectral density can be approximated by a single 
power law. Here we have taken into account that both 
magnetic and energy fluxes are almost constant along the 
magnetic field lines. In such a way the spectral param- 
eter T> on the solar surface can be evaluated by order 
of magnitude if we know the satellite spectral parameter 
p(sat) anc j .j.]^ ra y f constant components of the 
magnetic field 
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The observed power law for the energy density oc l// 2 
is theoretically explained in the framework of Langevin- 
Burgers MHD modeli 2 -. The ID calculations for the time 
and noise averaged spectral density of Alfven waves gen- 
erated by a white noise random driver for the external 
force density, modeling the influence of the convective 
stochasticity 

(F(t uZl )F(t 2 , z 2 )) = tp 2 5{t 1 - t 2 )S( Zl - z 2 ) (8) 

reveal the same inverse proportionality to the second 
power of the Alfven frequency 
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Comparison of the energy flux theoretically derived on 
the basis of the Langevin-Burgers approach applied for 
modeling the role of the turbulence for generation of 
Alfven waves with the experimentally observed energy 
flux can give us a reliable assessment for the Burgers pa- 
rameter r 
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Thus, if we consider Burgers approach as adequate for a 
turbulence model description, we can extract information 
for the turbulence spectrum in the photosphere, at the 
footpoints of the magnetic field lines 
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According to a recently proposed scenari o ' in the 
spirit of earlier ideas for wave heatin g 15 ' 16 Alfven waves 
serve as mediators, carriers of energy from the turbulent 
photosphere to the hot solar atmosphere, where in a small 
region the high-frequency waves attenuate intensively 
and heat the corona. RefJ^ treats the emission, while 
RefJ^ considers the bulk absorption of Alfven waves. In 
this work we present an evaluation for the absorbed in the 
transition zone energy flux, whose strong temperature de- 
pendence is determined by the temperature dependence 
of the cut-off frequency f c and naturally results in a sharp 
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The absorbed energy flux is taken from Voyager 1 mag- 
netometric data analysis, but it may also be obtained by 
a thermodynamical approach. If the comparatively small 
effects associated with radiative losses and compression 
are neglected, the absorption rate will be related only to 
the plasma internal energy density e and the conducted 
work 
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w ev + pv, e= -p, p = n e T e + n p T p , 

where p and v are respectively the plasma pressure and 
velocity. In such a way we can derive the approximate 
rate for the velocity of the solar wind, driven by sharp 
coronal temperature increase due to absorption of inten- 
sive high frequency Alfven waves, for which the transition 
region plasma is opaque 
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As the shear viscous friction heats the heavy 
particles , 10 ' 14 the proton temperature T p is significantly 
higher than that of the electrons T e , however, for an or- 
der of magnitude evaluation here we suppose the electron 
temperature T e to be similar to that of the protons T p . 



III. DISCUSSION AND CONCLUSIONS 

With a logarithmic accuracy we have derived an ex- 
plicit formula for the initial velocity of the solar wind. 
This formula Eq. (|14j) is completely based on experimen- 
tally accessible parameters and can be easily rejected if 
it gives more than 3 orders of magnitude difference. But, 
if this model remains valid, let us briefly discuss what 
has to be done as a future perspective in order to finally 
solve the perplexing longstanding mystery for the origin 
of coronal heating and solar wind. First of all, numerical 
simulations on MHD with a white noise random driver-S. 
(Langevin-Burgers MHD) have to be repeated to repro- 
duce a (5-like maximum of the energy dissipation density 
at the transition region. In other words all theoretical 
models qualitatively explaining the Voyager 1 data for 
the frequency dependent spectral density of the magnetic 
field have to be compared with random driver computer 
simulations of coronal heating. Confirmation of a narrow 
maximum of volume heating power due to self-induced 
plasma opacity is a routine task for further numerical in- 
vestigations. The numerical analysis could improve the 
present analytical evaluation incorporating, for example, 
the reflection of Alfven waves by the jump in plasma den- 
sity. 

Here we wish to insert a short historical remark. 
Stochastic mechanics in general was introduced by 
Langevir>l 7 in 1905 to explain the Brownian motion. 



Later, Burgers^ in 1948 introduced the white noise ran- 
dom driver in the hydrodynamics of turbulence. Much 
later, in 1995, PolyakoviS derived Kolmogorov power 
laws2°- using Langevin-Burgers approach, but his research 
remained unobserved in astrophysics, not to speak about 
heliophysics. That is why using random number genera- 
tor in MHD simulations is simply called random driver, 
but the magnitude of the noise is almost never evaluated 
by comparison to real experimental data. 

In the recent work we have only evaluated the area 
under the sharp maximum of dissipation and now it is 
time to perform state of the art calculations with a re- 
alistic random force, whose spectral density corresponds 
to the data taken by the satellites' magnetometers. For 
example, the turbulence spectrum can be treated with 
the based on the Langevin-Burgers model assessment 
Eq. (jlip. We propose that the physical mechanism for 
the origin of solar wind and coronal heating is already re- 
vealed and it is only a matter of honest numerical work 
to create a coherent picture. The results for the temper- 
ature jump and dissipation maximum have to be com- 
pared with the observations and other theoretical scenar- 
ios such as magnetic reconnectionsj^i for example. Each 
model unable to reproduce a sharp temperature increase 
in height has to be assigned to the waste basket. The 
same can be said for the models based on Ohmic heating 
which predict an electron temperature higher than that 
of the protons. Though the latter statement seems to 
be true for the X-ray bright points, it is definitely not 
valid for the corona as a whole. We are unaware, for 
example, how magnetic reconnection theory can explain 
why the proton temperature is higher than that of the 
electrons T p > T e . Analogously we have not found any 
reference explaining how turbulent cascade can lead to a 
very sharp increase of the coronal temperature. 

Second, in order to evaluate the properties of the solar 
surface as a random driver, all the old data from Voy- 
ager 1 has to be meticulously analyzed and a detailed 
investigation by the forthcoming Solar Orbiter mission 
has to be planned. 

Third, the Sun is a unique system for investigation 
of convective turbulence. It will be very interesting to 
compare the satellites' data for the magnetic field spec- 
tral density with results based on theoretical modeling 
of turbulence. The maturity of solar physics can stimu- 
late significant development of the achievements in con- 
temporary turbulence research. For instance, the l// 2 
power law by Burlaga and Mish£ corresponds to one- 
dimensional (k Xl k y — 0) propagation of Alfven waves in 
the framework of Langevin MHD,— as the whole noise is 
created by the random motion of the funnel foot-points. 

Next, absorbption of Alfven waves could be an im- 
portant mechanism in many cases of space plasmas in- 
cluding accretion disks as well. MHD waves in accretion 
disk amplification by shear flow is the main mechanism 
of transformation of gravitational energy into heat; for 
references see the preprint ^ The magnetized wind from 
an accretion disk should follow the magnetic force lines 
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and plasma concentrated above the disk center will create 
an accumulative jet perpendicular to the disk. In such a 
way disk jets could be recognized as a wind coming from 
magnetized turbulent plasma. 

Lastly, now we operate with a realistic 3D model for 
the distribution of the magnetic field from the solar 
surface to the satellite. The perturbation of magnetic 
field lines serves as a string of a harp to deliver the in- 
formation about the solar turbulence from the photo- 
sphere to the magnetometer. Owing to the propagation 
of Alfven waves we can "listen" to the sounds of the 
great solar symphony. Due to absorption of the high fre- 
quency modes, however, at the transition region where 
the temperature jump occurs we are able to hear only 
the basses, whereas for the ultra- violin band we remain 
absolutely deaf. A current problem which deserves to 
be put on the agenda is to observe the powerful high 
frequency Alfven modes under the coronal A kamikaze 
satellite could give some very important information, but 



for systematic research we need to learn how to extract 
the behavior of the solar surface as a random driver using 
optical data. We conclude that the first important step 
in this direction is to establish correspondence between 
satellite magnetometric data and Doppler shift spectra 
for some bright events on the solar surface. Only af- 
ter a proper incorporation of these ingredients we can 
conclude that our understanding of heating mechanism 
of the solar corona is complete and we have disclosed a 
very important case of heating of space plasmas. 
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